Astronomy <fc Astrophysics manuscript no. 5297'AA 


February 5, 2008 


(DOI: will be inserted by hand later) 





The X-ray spectrum of the bursting atoll source 4U 1728-34 

observed with INTEGRAL 

M. Falanga 1 ' 2 '*, D. Gotz 1 , P. Goldoni 1 ' 2 , R. Farinelli 3 , A. Goldwurm 1 ' 2 , S. Mereghetti 4 , A. Bazzano 5 , 

L. Stella 6 



O 
O 

(N 



CEA Saclay, DSM/DAPNIA /Service d'Astrophysique (CNRS FRE 2591), F-91191, Gif sur Yvette, France 
Unite mixte de recherche Astroparticule et Cosmologie, 11 place Berthelot, 75005 Paris, France 
Dipartimento di Fisica, Universita di Ferrara, via Paradiso 12, 44100 Ferrara, Italy 
INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica, Milano, via Bassini 15, 20133 Milano, Italy 
INAF-Istituto di Astrofisica Spaziale e Fisica cosmica, via del Fosso del Cavaliere 100, 00133 Roma, Italy 
INAF-Osservatorio Astronomico di Roma, via Frascati 33, 00040 Monteporzio Catone (Roma), Italy 



the date of receipt and acceptance should be inserted later 



> 

o 

m 
m 

o 

O 

6 

c3 



Abstract. We present for the first time a study of the 3-200 keV broad band spectra of the bursting atoll source 
4U 1728-34 (GX 354-0) along its hardness intensity diagram. The analysis was done using the INTEGRAL public 
and Galactic Center deep exposure data ranging from February 2003 to October 2004. The spectra are well 
described by a thermal Comptonization model with an electron temperature from 35 keV to 3 keV and Thomson 
optical depth, tt, from 0.5 to 5 in a slab geometry. The source undergoes a transition from an intermediate/hard 
to a soft state where the source luminosity increases from 2 to 12% of Eddington. We have also detected 36 type 
I X-ray bursts two of which show photospheric radius expansion. The energetic bursts with photospheric radius 
expansion occurred at an inferred low mass accretion rate per unit area of rh ~ 1.7 x 10 3 g cm -2 s _1 , while the 
others at a higher one between 2.4 x 10 3 - 9.4 x 10 3 g cm s . For 4U 1728-34 the bursts' total fluence, and the 
bursts' peak flux are anti-correlated with the mass accretion rate. The type I X-ray bursts involve pure helium 
burning either during the hard state, or during the soft state of the source. 
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1. Introduction 

> 

k> , Low-mass X-ray binary systems (LMXBs), consisting of 
■ an accreting neutron star (NS) and a main-sequence donor 

5^ ! s t ar (M< 1 M Q ), are usually classified as atoll or Z 
sources, according to the path they describe in an X-ray 
colour-colour diagram (CCD) or hardness-intensity dia- 
gram (HID). The main difference between atoll and Z 
sources is in the accretion rate, the magnetic field and 
probably in the spin rate and the geometry of the sys- 
tem. For atoll sources, the positions along the character- 
istic branches have historically been called the "island", 
"lower banana" and "upper banana" branches, which indi- 
cate the spectr al states probably as a function of the mass 
accretion rate l|Hasinger fc van der Klis. 1989). The mass 
accretion rate increases from the island state (hard spec- 
tral states) where the soft emission is much reduced and 
the spectrum is dominated at high energies up to 200 keV, 
to the banana states (soft spectral states) where most of 



the energy is emitted be low ~ 20 keV (e.g.. iBarreti l200ll 
lOierlihski fc DoneLl2002il . 

The atoll sources show luminosities smaller than those 
of Z sources (generally less than 10 per cent of Eddington, 
-^Edd)) low magnetic fields (B^ 10 8 G), and spectral vari- 
ability, or equivalcntly, motion along the banana branch 
over a timescale of hours to days, and more slowly in the 
island state from days to weeks. 

Some atoll sources also show time variability prop- 
erties in their power spectra, such as Quasi-Periodic 
Oscillations (QPOs), which are strongly correlated with 
the p osition on the CCD/HID l|Hasinger fc van der Klisl. 
1989). Almost all the timing features characterizing 
the power density spectra of these sources seem to 
vary in a smooth and monot onic way when the source 
move s along its CCD/HID llWiinands fc van der Klisl 
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Il999t IPsaltis. Belloni fc van der Klisl Il999j) . implying a 
tight correlatio n between spectral a nd temporal behaviour 
(for a review see lvan der Klisl 12004]) . In most atoll sources, 
QPOs are observed in the frequency range from millihertz 
to the kilohertz. Kilohertz QPOs or twin kHz QPOs in the 
frequency range from 300-1200 Hz are usually observed in 
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the i sland state and lower banana branch l|van der Klisl 
2004 and references therein). 

In some of the ato ll sources, type I X -ray bursts 
are also observed (e.g.. iLewin fe Josa Il983[) . Type I X- 
ray bursts are characterized by a fast rise and expo- 
nential decay with durations ranging from seconds to 
tens of minutes. These bursts are due to unstable hy- 
drogen/helium burning in a thin shell on the NS sur- 
face ( see e.g., the review bv lLewin. van Paradiis. fc Taaml 
Il993l) . The burst spectrum can be described by black- 
body radiation with cooling during the decay of the 
burst. The gradual decay persists longer at lower en- 
ergies, indicating that the burst spectrum is charac- 
terize d by steadily decre asing blackbody temperatures 
(e.g., ILewin fc Josl Il983j) . During these bursts nearly- 
coherent oscillations are sometimes observed, the fre- 
quencies of which are in the rather narrow range be- 
tween 300 and 600 Hz JStrohmaver fc Markwardtl ll 



19811 iBasinska et al l Il984t IWhite et all Il98fit 
Claret et all Il994t ISchulzl. Il999^ and^ more re- 



Tz llStrohm 

IStrohmaver fc BildstenL |2004). This frequency is inter- 
preted as the NS rotation frequency due to a hot spot 
(or spots) in an atmosph eric layer of the rotating NS 
1 Chakrabartv et all Eool . 

The CCDs diagrams and HIDs are powerful ways of 
parameterizing the spectral changes using physically moti- 
vated spectral models to understand the underlying phys- 
ical changes in the source emission. In this work, we study 
for the first time the broad band (3-200 keV) spectral 
shape of the atoll source 4U 1728-34 as a function of its 
HID position. The correlation of type I X-ray burst prop- 
erties with source state in the HID is also investigated. 
We analyze the International Gamma- Ray Astrophysics 
Laboratory {INTEGRAL) 2003/2005 data concentrating 
also on the high energy emission of 4U 1728-34. This pa- 
per deals only with type I X-ray bursts, therefore whenever 
we write burst (s), we mean type I X-ray burst (s). 

1.1. The source 4U 1728-34 



The persistent LMXB source 4U 1728-34 (or GX 354-0) 
is a well known burster discovered by Uhuru and classi- 
fied as an atoll sour ce from its colour-colour diagram and 
timing proper ti es llForm an., Tananbaum Jonesl [l 9 76 ; 
iHoffman et all Il97fil lHasinger fc van der Klisl FoSfllv . 
However, little is known about the system as the opti- 
cal counterpart is not identified due to the high optical 
extinction in the direction of the Galactic center. The dis- 
tance is poorly known, the current estimate being between 
4.4-5.1 kpc using photospheric r adius expansion burst 
luminosity as a stan dard candle (|Di Salvo et all [2000: 
iGallowav et"ail l2003|) . This is based on the assumption 
that the bolometric burst peak luminosity during photo- 
spheric radius expansion is saturated at the Eddington 
limit. 

The time averaged persistent soft (0.1 keV) to 
hard (300 keV) X-ray emission spectrum of 4U 1728- 
34 was studied in the past using Einstein, SAS 3, 
EXOSAT SIGMA and ROSAT data llGrindlav fc Hertzl 



cently, using R X TE, B e ppoSAX, ASCA and Chandra 
(Di S alvo et al 1 120001 Piraino. Santangclo & Kaaret, 
| 2000l iNarita. Grindlav fc Barretl l200ll iD'Ai et al. . 
2006J). The spectrum was well fitted by a thermal 
bremsstrahlung or thermal Comptonization model as 
well as a Gaussian emission line at the energy ~ 6.7 
keV. Different spectral parameters were found for the 
electron temperature and Thomson optical depth, which 
can be explained by the varying state of the source in its 
CCD/HID position. 

During all spectral stat es, bursts of ~ 10-20 s du- 
ration were observed (e.g., ICornelisse et all l2003h . The 
burst spectrum can be fitted by a blackbody with tem- 
perature increasing from 1 to 3 keV during the burst rise 
time. Most of these bursts show evid e nce fo r photospheric 
radiu s expansion ( Hoffma n et al 1 Il976t IWhite et al " 



1983 iDav fc Tawa.raUl900l ^Foster. Ross, fc Fabian! ll 98f : 
Franco, 200l| Ivan Straatenl Mm ICallowav et all |2003, 



Flux oscillations at 363 Hz during the bursts 
were discovered in this source with the Rossi X- 
ray Timing Explorer (RXTE), and interpreted as the 
spin frequency o f the NS l|Strohmaver et all Il996t 
IChakrabartv et all l2003h . Apart from the bursting ac- 
tivity, 4U 1728-34 also exhibits complex behaviour on 
short time scale. It was also the first neutron star LMXB 
to show twin kHz QPOs in the persistent emission 
with a fr equency separation alm o st constant at ~ 350 
Hz ( e .g.. iMendez fc van der Klisl Il999| iDi Salvo et all 
l200ll iMigliari. van der Klis. fc Fender! 120031 and refer- 
ences therein). 



2. Observations and Data 

The present data set was obtained with the INTEGRAL 
observatory l|Winkler et all 120031) using the 2003/2005 
Core Program Galactic Center Deep Exposure observa- 
tions and publicly available data. We analyzed data from 
the Imager on Board the Integral Satellite (IBIS) coded 

mask INTEGRAL Soft Gamma-Ray Im ager (ISGRI) 
• 1 1 73 I 1 r"^ — ~i 

(Ubertini ct al., 2003; Lcbrun et al., 2003) at energies be- 
tween 20 an d 200 keV and from the JEM-X monitor, mod- 
ule 1 and 2 l|Lund et alll2003|) between 3 and 20 keV. For 
ISGRI and JEM-X, the data were extracted for all point- 
ings within < 3?5 from the source direction. The available 
data consist of 356 individual pointings for a total effec- 
tive exposure of 640 ks. Data reduction was performed 
using the standard Offline Science Analysis (OSA) soft- 
ware version 5.0 dist ributed by the INTEGR AL Science 
Data Centre (ISDC) l|Courvoisier et allbOO^ . The algo- 
rithms used in the IBI S /ISGRI analysis are described in 
lOoldwurm et~aTI l|2003l) . 

The brightest bursts were detected using the 
INTEGRAL Burst A lert System (IBAS) software 
llMereghetti et all l2003|) running at ISDC. IBAS is ded- 
icated to the real time discovery and localization of 
gamma-ray bursts, transient X-ray sources and bursts 
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in the IBIS/ISGRI data stream. Bursts are identified by 
scanning the data binned using different timescales (from 
a few milliseconds to 100 s) and energy ranges. One of 
IBAS operation modes, running on the 15-40 keV energy 
interval and 10 s time scale, is particularly suited to de- 
tect type I X-ray bursts. Typically the bursts are localized 
with a ~ 3' uncertainty. The fainter bursts were detected 
using 2 s rebinned JEM-X light curves in the 3-20 keV en- 
ergy band. Most of the bursts detected with JEM-X were 
also detected with ISGRI in the 18-40 keV energy band, 
these flux excesses were confirmed also by the imaging as 
related to the source position. So, we call the bursts de- 
tected with IBAS "bright bursts" and the ones detected 
only with JEM-X "weak bursts" . 

The burst light curves are based on events selected ac- 
cording to the detector illumination pattern for 4U 1728- 
34; for ISGRI we used an illumination threshold of 0.6 for 
the energy range 18-40 keV. For JEM-X we used the event 
list of the whole detector in the 3-6 keV, 6-12 keV and 
12-20 keV energy band. In Table ^ we report the burst 
start time detected with INTEGRAL from February 2003 
to January 2005. The start time for each burst was de- 
termined when the intensity was 10% of the peak above 
the persistent intensity level. Note that six of the IBAS de- 
tected bursts were outside of the JEM-X field of view. Due 
to the lack of the low energy data we excluded these bursts 
from the analysis. However we report the burst start times 
and durations (see Table QJ. 

3. Results 

To study the light curves and spectra of 4U 1728-34 we 
first deconvolved and analyzed separately the 356 point- 
ings and then combined them into a total mosaic image 
in the 20-40 keV energy band. In the mosaic 4U 1728-34 
is clearly detected, as well as the nearby source the Rapid 
Burster (4U 1730-335). The source position offsets with 
respect to the catalog positions are 0.'2 for 4U 1728-34 
and 0.'23 for the Rapid Burster. This is within the 90% 
con fidence leve l assum ing the source location error given 
bv iGros et al.1 l)2003|) . The derived angular distance be- 
tween the two sources is ~ 31. '5. Due to the fact that 
INTEGRAL is able to image the sky at high angular res- 
olution (12' for ISGRI and 3' for JEM-X), we were able 
to clearly distinguish and isolate the high-energy fluxes 
from the two sources. This allowed us to study the X- 
ray emission of 4U 1728-34 during its entire observation 
without contaminat ion (see also for the Rapid Burster 
iFalanga et all l2004|) . In Fig. [2 we show the JEM-X (top) 
3-20 keV and ISGRI (bottom) 20-100 keV light curves 
extracted from all the 356 images. Pointings during which 
bursts were detected are marked in red (bright burst) and 
green (weak burst). The source is very active and strong 
variations on a short time scale are not uncommon. Bursts 
were observed during all intensity states. We also observed 
an isolated strong flare (around 53058 MJD), during which 
the source count rate increases from ~ 5 cts s _1 to ~ 50 
cts s _1 in the soft energy band, 3-20 keV, and shows a 
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Time since MJD 52698 (days) 

Fig. 1. JEM-X and IBIS/ISGRI light curves of 4U 1728- 
34 in the 3-20 keV and 20-100 keV energy bands. We 
indicate with red and green the net intensity (burst sub- 
tracted) emission where bursts were emitted. Each point 
corresponds to a pointing of ~ 1800 s duration. The 21- 
day flare is shown on a smaller time scale in the insets. 

different trend at higher energy, 20-100 keV. During this 
flare, lasting ~ 21 days, the source undergoes a transition 
from a hard to a soft state (see Sect. 13. 3|) . Note that dur- 
ing this flare two weak bursts happen at the peak of the 
3-20 keV flux. 

3.1. Burst light curves 

In total, we observed 36 bursts, which we label numeri- 
cally according to time. Their main properties are listed 
in Table □ Based on JEM-X/ISGRI light curves in the 
3-6 keV, 6-12 keV, 12-20 keV and 20-40 keV energy 
bands, rebinned at 0.5 s, we determined that the bursts 
have rise times of 1 ± 0.5 s. The decay time at 12-20 keV 
energy is significantly shorter than that in the 3-6 keV 
range, indicating spectral softening during each burst de- 
cay. This was also verified by determining hardness ra- 
tios, using the three JEM-X energy bands. For JEM-X 
the flux decays to quiescence with an mean exponential 
time-scale F oc e - */ 4 - 5 (see Table H). At high energy, 
the bursts were only significantly detected in the first 
AT Bu rst ~ 5-10 s (see Table [lj. In Fig. |21 we show six 
JEM-X/ISGRI representative burst light curves in dif- 
ferent energy bands. For bursts 2 and 4, a double peak 
profile is evident at high energy (lower panel) within the 
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Table 1. Log of INTEGRAL burst observations analyzed in this paper and spectral properties of the bursts. 



Burst 


IBIS/JEMX 


IBIS/JEM-X 


IBIS 


JEMX 




Burst 


Parameters 






Nr. 


Date 


T start (UTC) 


ATflurst 


» exp 


fn 


E b 


Fpcak 


T 


7 


a 


it 


2003 Feb 28 


07:55:01 


8.4(3) 
















2+ 


2003 Mar 01 


00:04:45 


8.2(3) 


6.8(5) 


1.7(2) 


8.3(3) 


9.6(3) 


8.6(6) 


0.012(4) 


~ 110 


3+ 


2003 Mar 01 


16:05:27 


9.4(3) 


_ 




_ 


_ 





_ 


_ 


4t 


2003 Mar 02 


07:42:17 


11.0(3) 


6.8(5) 


1.7(2) 


8.9(3) 


9.9(3) 


9.0(4) 


0.016(3) 


~ 100 


5 


2003 Mar 21 


06:09:48 


5.7(3) 


4.5(5) 


3.2(2) 


3.5(3) 


3.9(3) 


9.0(9) 


0.07(7) 


_ 


6 


2003 Mar 21 


15:46:38 


3.5(3) 


3.6(5) 


4.3(2) 


3.2(2) 


7.9(2) 


4.1(3) 


0.051(2) 


~ 270 


r 


2003 Apr 03 


08:40:16 


4.2(3) 
















8+ 


2003 Apr 06 


19:45:28 


5.8(3) 


5.9(5) 


3.7(2) 


5.1(4) 


7.4(4) 


6.9(6) 


0.046(8) 




9 


2003 Apr 12 


07:12:20 




3.2(5) 


3.5(2) 


2.6(2) 


4.4(2) 


5.8(5) 


0.073(6) 




10 


2003 Apr 14 


09:05:15 




3.4(5) 


3.7(2) 


1.9(2) 


3.4(2) 


5.5(7) 


0.103(8) 




11 


2003 Apr 14 


20:13:47 




5.0(5) 


4.5(3) 


2.7(3) 


3.8(2) 


7.1(9) 


0.111(8) 


~ 610 


12 


2003 Apr 21 


09:31:41 




5.4(5) 


3.7(2) 


2.1(2) 


3.6(2) 


5.8(6) 


0.062(7) 




13 


2003 Apr 19 


22:57:35 


3.5(3) 


3.4(5) 


2.7(2) 


3.1(2) 


5.1(2) 


6.1(5) 


0.049(4) 




14 


2003 Aug 23 


16:14:00 


4.2(3) 


5.0(5) 


2.7(2) 


3.5(2) 


6.1(2) 


5.7(4) 


0.042(4) 




15 


2003 Aug 25 


18:24:08 


6.0(3) 


4.1(5) 


3.9(2) 


3.6(2) 


6.6(2) 


5.5(3) 


0.056(3) 




16 


2003 Aug 28 


01:23:58 


7.6(3) 


3.6(5) 


3.2(5) 


3.9(5) 


5.5(5) 


7.1(6) 


0.055(5) 




17 


2003 Aug 28 


06:01:29 


4.6(3) 


2.6(5) 


3.3(5) 


3.1(5) 


5.3(5) 


5.8(6) 


0.057(6) 


~ 160 


18 


2003 Aug 31 


15:54:14 


6.8(3) 


6.5(5) 


3.1(4) 


3.1(4) 


6.3(4) 


5.0(7) 


0.046(7) 




19 


2003 Aug 31 


20:00:09 




4.0(3) 


3.5(3) 


1.9(4) 


3.9(4) 


4.9(9) 


0.08(1) 


~ 240 


20 


2003 Aug 31 


23:48:20 




5.6(3) 


4.6(2) 


3.1(3) 


3.7(3) 


8.4(9) 


0.12(1) 


~ 200 


21 


2003 Sep 03 


03:26:33 


5.0(3) 


3.1(5) 


3.1(3) 


3.8(5) 


7.1(5) 


5.3(8) 


0.04(7) 


_ 


22 


2003 Sep 03 


08:39:29 


6.8(3) 


4.9(5) 


2.9(5) 


2.8(3) 


4.5(3) 


6.2(7) 


0.06(5) 


~ 170 


23 


2003 Sep 03 


13:07:21 


2.7(3) 


3.9(5) 


2.8(4) 


2.1(4) 


4.4(4) 


4.8(9) 


0.06(1) 


~ 220 


24 + 


2003 Oct 08 


09:58:35 


5.5(3) 


_ 


_ 


_ 


_ 


_ 


_ 


_ 


25 + 


2004 Feb 17 


04:47:27 


9.9(3) 


_ 


_ 


_ 


_ 


_ 


_ 


_ 


26+ 


2004 Feb 19 


21:06:41 


5.6(3) 


1.97(5) 


2.4(4) 


4.2(4) 


8.2(4) 


5.1(5) 


0.028(5) 


_ 


27 + 


2004 Feb 20 


12:00:21 


6.3(5) 


_ 


_ 


_ 


— 


— 


_ 


_ 


28 


2004 Mar 08 


04:14:44 


4.5(3) 


5.2(5) 


12.4(3) 


3.3(3) 


6.8(3) 


4.8(6) 


0.171(5) 


— 


29 


2004 Mar 09 


00:58:20 


4.9(3) 


3.6(5) 


10.4(4) 


3.5(4) 


6.9(4) 


5.0(7) 


0.14(5) 


~ 690 


30 


2004 Sep 07 


21:30:45 


5.0(3) 


3.1(5) 


3.8(3) 


4.2(3) 


7.0(3) 


6.0(6) 


0.056(5) 


- 


31+ 


2004 Sep 08 


12:44:37 


6.8(3) 


4.9(5) 


4.1(3) 


4.6(3) 


7.3(3) 


6.3(6) 


0.052(5) 


~ 450 


32 f 


2004 Oct 02 


05:59:14 


2.7(3) 


3.9(5) 


4.1(3) 


4.1(3) 


6.9(3) 


5.9(6) 


0.054(6) 




33 f 


2004 Oct 02 


21:51:07 


6.1(3) 


5.5(5) 


3.7(4) 


3.6(4) 


7.7(4) 


4.7(5) 


0.074(5) 


~ 550 


34 + 


2004 Oct 04 


12:29:37 


4.5(3) 


5.2(5) 


4.7(3) 


3.9(3) 


7.3(3) 


5.3(5) 


0.059(6) 




35+ 


2004 Oct 05 


05:59:09 


6.0(3) 


4.6(5) 


5.5(3) 


4.0(3) 


6.7(3) 


5.9(5) 


0.082(5) 


~ 780 


36 f 


2004 Oct 11 


12:59:09 


4.9(3) 


6.8(5) 


1.7(3) 


3.5(3) 


7.0(3) 


5.0(5) 


0.091(6) 





Notes : ' bright type I X-ray bursts. * not a statistically significant detection in ISGRI. The burst exponential decay time, r oxp , 
is in units of seconds in the 3-20 keV energy band. The mass accretion rate per unit area of the persistent emission before the 
bursts, rh, are expressed in units of 10 3 g cm -2 s _1 , assuming a canonical NS radius of 10 km. The net burst fluence, Eh, and the 
net peak flux, -Fp ea k, (both persistent emission subtracted) are in units of 10 -7 erg cm -2 and 10 -8 erg cm -2 s _1 , respectively. 
The bursts parameters r = i?b/-£pcak is in units of seconds, and 7 = FpcrsZ-Fpcak is also listed for each burst. The quantity a 
is defined as the ratio of the total energy emitted in the persistent flux to that emitted in bursts, a — F pcls At/ Eh, where At 
is the time interval from the beginning of one burst to the next. We report only the most confident a values, i.e. without data 
gap between two burst. 



first 10 s, while during this time the intensity at lower 
energy (upper panel) remains constant. This can be inter- 
preted as a consequence of a photospheric radius expan- 
sion (P RE) episode during the first part of the outburst 
(see e.g. JVerbunt fc van den Heuveill99RT) . When a burst 
undergoes a PRE episode, due to radiation pressure, the 
luminosity, L, remains constant at the Eddington value, 
^Edd, and, at high energy, a double peak profile can be ob- 
served. Such behaviour is well known for 4U 1728-34 and 
was observed previously in different bursts and recently 



reported with RXTE data llFrancol l200li | van Stra atcn. 
l200lt iGallowav et all EofH l200rt . 



We attempted to test the hypothesis of a PRE us- 
ing the four different JEM-X/ISGRI energy bands and we 
measured such a double peak behaviour in the hardness 
ratio only for burst 2 and 4. 
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Fig. 2. Some examples of bright type I X-ray bursts detected from 4U 1728-34. For each burst the JEM-X (3-20 keV, 
upper panels) and IBIS/ISGRI (18-40 keV, lower panels) net light curves are shown (background subtracted). The 
time bin is 0.5 s for both IBIS/ISGRI and JEM-X. Bursts 2 and 4 show evidence of photospheric radius expansion. 



3.2. Hardness Intensity Diagram of the persistent 



emission 



To generate a CCD (see for this sourc e e.g., 
Hasinger fe van der Klisl Il989t iDi Salvo etail l200ll 
Francol 1200*1^ we used the background subtracted and 
burst subtracted JEM-X light curve with a 500 s time 
resolution. We defined the soft colour as the logarithm 
of the ratio of the count rate in the energy range 3-4.8 
keV to the count rate in the range 4.8-7.7 keV. Similarly, 
we defined the hard colour in the range 7.7-12.5 keV and 
12.5-20 keV. In this case the CCD was inadequate as 
the colour errors were as large as the colour variation. In 
order to reduce the errors on the colours, we generated 
a HID based on the net count rates in the 3-6 keV and 
6-20 keV obtained in the single individual pointings of 
~ 1800 s each. In our case the source behaviour is much 
better traced in the HID rather than in the CCD, because 
the statistical uncertainties are reduced along one of the 
two axes. The HID is shown in Fig. |3 with the positions 
superposed where the bursts were detected. Note that the 
patterns observed in a CCD are often also recognizable in 
the corresponding HID obtained by replacing soft colour 
by intensity, as applied e.g. on the Z sources GX 17+2 



Pi Salvo et all 1200(1 iFarinelli et all l2005h . GX 5 



Muno. Rcmillarc 



Lie 

ard fc Chakrabartl l2002t) and on the ato ll 
1 l|Muno. Remillard fc Chakrabartl 12002ft . 



-1 



source Aql X- 
This results from the intensity being dominated in the 
lower bands by soft photons whose thermal component 
correlates with temperature and therefore colour. The 
HID could indicate the banana branch of the atoll, with 
inferred accretion rate increasing along the atoll from 
left to right on the upper branch to the lower banana 
branch. To study the spectral variability of the source 
as a function of the HID position we divided the HID 
into nine regions, as shown in Fig. The intervals were 
chosen in order to have good statistics for the spectral 
analysis. We averaged the different intensities with nearly 
the same hardness ratios and considering the nearness 
colours these spectra should be very similar. We also 
studied in more detail the 21 day long flare (see Fig. HJ, 
divided among the boxes 6-9, represented by the filled 
dotted points in Fig. The HID of the flare shows a 
transition from a hard to a soft state during 21 days with 
an intensity ratio / m ax/-^min ~ 10 (see Sect.0J. This trend 
held also for the persistent emission, which also moved 
from a hard to soft state (boxes 1-5). This overall trend 
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3.3.1. Persistent emission variability 
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Fig. 3. HID of 4U 1728-34. The hardness is the ratio of the 
count rates in the 6-20 keV to 3-6 keV and the intensity 
is the 3-20 keV count rate. Each point (burst-subtracted) 
corresponds to ~ 1800 s of integration time. A typical 
error bar is shown at the top left corner. The boxes rep- 
resent our selection of the data for the spectral analysis. 
The choice of the boxes (1-5 and 6-9) was made in or- 
der to study the spectral evolution along the atoll and 
to compare our result with previously studied 4U 1728-34 
spectra. The filled dotted points are from the flare with 
-fmax/^min ~ 10 in the 3-20 keV energy band which lasted 
21 days (see Fig.JTJ). The bursts detected with JEM-X and 
ISGRI are also indicated. Bursts 2 and 4 show evidence of 
a photospheric radius expansion episode. 



can also be seen in the RXTE/ ASM light curve (Fig. |7J) 
for the period of this observation. 



We studied the spectra of 4U 1728-34 along the atoll pat- 
tern for the different boxes reported in Fig. [3] A broad 
band 3-200 keV spectrum was obtained for each box us- 
ing the joint JEM-X/ISGRI data and removing the time 
intervals corresponding to the bursts. 

Atoll source spectra are successfully fitted by a two- 
component model consisting of a multicolour blackbody 
soft X-ray emission and a Comptonized spectrum, for 
the hard X-ray emission. The emission from the accre- 
tion shock on a neutron star is expected to be produced 
by thermal Comptonization of soft seed photons from the 
star, or from an optically thick b oundary layer between 
accretion disc and NS surfac e (e.g.. lZeldovich fc Shakural 
Il969t lAlme fc WilsonL Il973l) . We model the sh ock emis- 
sion by the COMPTT model ijTitarchukl Il994h which is 
an analytic model describing Comptonization of soft pho- 
tons in a hot plasma. We approximate the accretion shock 
geometry by a plane-parallel slab at the neutron star sur- 
face. The main model parameters are the Thomson optical 
depth across the slab, Tt, the electron temperature, kT e , 
and the seed photon temperature, fcr soc( j. The soft thermal 
emission, kT so f t , believed to be produced b y the disc, is fit- 
ted b y a multi-temperature disc model of ijMitsuda et all 
1984). We were not able to constrain the interstellar col- 



umn density ATjj (as the JEM-X bandpass starts above 
3 keV), so we have set iVn to the value 2.5 xl0 22 cm~ 2 
observed by ROSAT and BeppoSAX. We add to the fit 
a Gaussian line around 6.4 keV with a fixed width at 
a = 0.34 keV as observed previously (e.g.. iDi Salvo et all 
2000t lPiraino. Santangelo fc Kaareti 120001) . but it was not 
statistically required. The multi-temperature disc model 
was also not required in the fit during the low/hard state. 
Soft emission (presumably from the accretion disc) during 
this state is difficult to detect since the JEM-X bandpass 
begins at 3 keV. Compton reflection between 10-30 keV 
was not required by the fit. The best fit parameters for 
each box are reported in Tableland in Fig. [3] For all the 
fits, the normalization constants of the JEM-X response 
were within 1.08 ±0.09. The nine unabsorbed vF v spectra 
and the residuals of the data to the model are shown in 
Fig.H 



3.3. Spectral Analysis 

The spec t ral an alysis was done using XSPEC version 11.3 
l Arnaudl ll996T) . combining the 20-200 keV ISGRI data 
with the simultaneous 5-20 keV JEM-X data. A constant 
factor was included in the fit to take into account the un- 
certainty in the cross-calibration of the instruments. The 
factor was fixed at 1 for the ISGRI data. A systematic er- 
ror of 2% was applied to JEM-X/ISGRI spectra which 
corresponds to the current uncertainty in the response 
matrix. All uncertainties in the spectral parameters are 
given at a 90% confidence level for single parameters. For 
the s ource distance we use 4.5 kpc 1 Galloway et all fepOSl 
2f 



3.3.2. X-ray burst spectra 

To study the net bursts' spectra, we first extracted the 
persistent emission spectrum excluding the outburst in- 
terval. We verified that during each burst pointing the 
count rate was stable, then we used the same fit mod- 
els as described in Sect. 13.3.11 to fit the data. In Table 
^ we report the mass accretion rate per unit area, to, of 
the persistent emission for each burst pointing. This is 
given by to = M/A acc , where A acc is the area covered 
by 4irRfj S . We assume a 10 km radius for the neutron 
star. The persistent mass accretion rate, M, is calculated 



from M 



7 pcrs 



r) c , where 77 <~ 0.2 is the accretion 
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Fig. 4. INTEGRAL/ JEM-X (3-20 keV) and ISGRI (20-200 keV) unfolded spectra of 4U 1728-34 along the "banana" 
track with the best fit diskbb and COMPTT models. Residuals between the data and model are shown in the bottom 
panel in units of sigma. The spectra are labeled corresponding to each box number selected in the HID in Fig. 03 The 
best fit parameters are shown in Fig. [5] and Table [21 



efficiency for a NS. The bursts occurred at a persistent 
emission level L pers ~ 2-8% I/Edd- 

We then studied the net bursts' time averaged spec- 
tra using for each burst the persistent emission spectra as 
background. The low statistics of JEM-X/ISGRI do not 
allow us to study the spectral evolution during the bursts. 
The spectra were extracted over the whole duration of 
each burst as determined from the JEM-X light curves. 
We also included in the fit the ISGRI spectra in the 18- 
40 keV energy bands. All obtained spectra were fitted by 
a photoelectrically absorbed blackbody. The unabsorbed 
flux between 0.1-200 keV, Fbd.bb, was calculated accord- 



ing to F\ 



bol.bb 



= 1.0763 x lO- 11 ^ 4 ^ erg cm- 2 s 



where Tbb is the blackbody temperature in unit of keV 
and Kbb(D,RNs) is the normalization of the blackbody 
component as returned by the fitting program. The black- 



body temperature of the time averaged burst spectrum, 
kTbb, was between 1.8 keV and 2.4 keV. The burst flu- 
ence, Eh, is calculated by integrating the measured Fboi.bb 
over the burst duration. The burst peak flux, F'peak, was 
calculated by converting the peak count rate to flux using 
the averaged spectrum of each burst. Taking into account 
the typ ical blackbody temp erature variation within the 
bursts l|van Straaterl l200l[) . the associated error on our 
conversion factor is of the order of 12%, which is within 
our statistical errors. The best fit results and calculated 
burst parameters are reported in Tabled 
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Table 2. Fit parameters of the persistent spectra in the energy range 3-200 keV 
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Iv soft 

kT seed (keV) l-18±g;S 

kT c (keV) 34.9lg 4 
~ n c;+0.l6 
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Fig. 5. Spectral parameters from Table [21 as a function of 
the position on the HID. 

4. Discussion 

4.1. Origin of X-ray emission 

We have fitted the 3-200 keV spectra of the atoll source 
4U 1728-34 along its track on the HID by a two-component 
model: a thermal Comptonization model together with 
a soft component. The spectra and the parameters are 
shown in Fig. ^ and Fig. |SJ respectively (see Table [2J . 
The change in the spectral shape between the hard state 
(box/spectrum 1) and the soft state (box/spectrum 5) is 
clearly visible, as also during the flare (box/spectra 6-9), 
probably caused by a dramatic change in the accretion 
flow geometry. During the spectral evolution, spectra 1— 

5, at low intensity the optical depth was tt ~ 0.5 with 
a high plasma temperature kT c ~ 35 keV (box 1), and 
at higher intensity the optical depth was tt ~ 1.4 with a 



plasma temperature as low as kT c ~ 6 keV (box 5). During 
this spectral change the luminosity increases from ~ 0.02 
to 0.05 LEddj where the seed photons temperature was al- 
most constant at kT S00( i ~ 1.3 keV. The same behaviour is 
observed also during the flare event (box 6-9) where the 
plasma temperature decreases from ~ 18 to 3 keV and 
the optical depth increases from 1 to 5 and the spectral 
change the luminosity increases from ~ 0.06 to 0.12 I/Edd- 
We observe a typical atoll spectral change from a low or 
intermediate/hard state to a high/soft state, where the 
accretion rate increases from the hard state to the soft 
state. The soft component could be associated with the 
radiation from the accretion disc. 

With our coverage of luminosity variation, 
F max /F m i n ks 2.5, observed during ~ 7.4 days inte- 
gration time between 2003 February and 2004 October 
we are probably not observing the full atoll track, even 
during the flare (see Fig. ■ We observe the source most 
likely from a intermediate/hard state, upper banana, 
to the soft state, lower banana state, and probably 
not during the island state. If the source was in the 
island state, we could also expect that more then two 
bursts with PRE are observed (see Sect. I4.2TI . The hard 
component observed during the hard/intermediate state 
is explained with Comptonization of seed photons from 
a NS surface by hot electrons in an inner optically thin 
accretion flow (or outer boundary layer). The optically 
thin Comptonization medium has temperature of > 20 
keV. The observed plasma temperature and optical depth 
during the spectra 1-3 are typically hard/intermediate 
state parameters, only during the flare we observe 
typically lower banana state parameters. In the lower 
banana branch spectra the Comptonized component is 
much softer, with temperature of ~ 6 keV and optical 
depth of 3> 1 as observed during the spectra 5-9. 

From our fits it results that the hard spectra can be de- 
scribed by unsaturated Comptonization of soft photons in 
the hot kT c - 20 - 30 keV optically thin r - 0.5-1 
plasma. However, different models have been proposed 
to explain the different emission region and emission 
mechanism during the different spectral states and evo- 
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lution llMitsuda, et all Il989t IZeldovich fc Shakvlral Il969t 
lAlme fc Wilsod.ll973tlw"hite. Stella fc Parmeii Il98^ ~ 

If the soft photons are emitted from a disc, then the 
inner disc radius can be estimated from the diskbb model 



as 



Rm 



2.7 4.5 



/col 

~L8 



0.5 



COS I 



1/2 



km, 



(1) 



where D is the source distance, i is the inclination an- 
gle of the system, f nr ,i is the ratio of the colo ur to ef- 
fective temperature l|Shimura fc Takaharal fl995) and rj is 
the correc tion factor for the inn er torque-free boundary 
condition I Gierlihski et all Il999l 77 = 2.7 for R ul ~ 7?ns 
and less for higher R[ n ). With D = 4.5 kpc, / co i = 
1.8 and using a mean inclination angle of i = 60° 
we find i? in = 15 ± 3 52 ± 13 km. For 4U 1728-34 
the observ ed QPO are related on their position in the 
CCD /HID llMendez fc van der KlisUl999tlDi Salvo et all 
l200lUMigliari. van der Klis. fc Fended. l2003|) . Among the 
most frequently used QPO models, the upper kHz QPO 
is related to the Kepler frequency, i.e. ^qpo = t'Kopicr ~ 
1184 Hz (i? in /15 km)- 3 / 2 Mg^ 4 , at the inn ermost stable 
orbit of the accretion disc (for a review see Ivan der Klisl 
120041 and reference therein). Using our measured inner 
disc radius, i?i„, the upper kHz QPO has to be in the fre- 
quency range from 1160 ± 300 Hz to 213 ± 53 Hz. These 
values are c onsistent with the up per kHz QPO values of 
4U 1728-34 (|Di Salvo et alll200l|) . However, using this re- 
lation we are not able to exactly locate the source position 
in the CCD/HID. 

We have measured the plasma temperature, op- 
tical depth, and flux of 4U 1728-34 over a range of 
positions on the HID for the first time. We com- 
pare with previously r eporte d time averaged spec- 
^^^^^^^^^^^^^r^^^^^^^n ^^^^t^ ^^^^^^^^^^^^^^^^^ ^^^^n 



White et all 


198a IClaret et all Il994l ISchulzL Il99£ 


Di Salvo et al. 


120001 IPiraino, Santangelo & Kaaret 



al 



2006). Accounting for the fact that these authors have 
reported time averaged spectra sometimes over possibly 
different spectral states, and sometimes only in energy 
ranges less than 20 keV, and insofar as the models used to 
fit the data are consistent, we found that all the previous 
reported spectral parameter s are consist e nt w ith our 
measurements. For example, IClaret et alJ l|l994|) . using 
SIGMA (30-200 keV) data, measured plasma tempera- 
tures KT C of ~ 5 and ~ 30 keV in two different source 
states. These correspond to our reported soft and hard 
states. Also using a broad band spectrum, BeppoSAX 
observed the source with kT e ~ 7 keV and t ~ 5, also 
consistent with an intermediate soft state. We would 
like to point out, in order to determine the source states 
with more confidence, how important it is for 4U 1728-34 
to have: (i) a high energy spectrum above 20 keV, and 
(ii) not to average the spectra over a long time period. 
(i) using the JEM-X and ISGRI broad band spectrum 
(box 1) from 3-200 keV, in the intermediate/hard state, 



the thermal Comptonization plasma temperature was ~ 
35 keV, considerably higher than the value (~ 10 keV) 
derived from the low energy 3-20 keV JEM-X data alone. 
(ii) a flux variation of a factor 1.12 will already change 
the spectral parameters and a factor 2 will drastically 
change the source state (see Table Note also, e.g., 
that during the flare we have I ma , x /I m in ~ 10 in the 
3-20 keV band which corresponds to a flux variation of 
-Fmax/-Fmin ~ 2 in the 3-200 keV band. This shows again 
for 4U 1728-34 how important the energy contribution is 
above 20 keV. 



4.2. Type I X-ray bursts 




t PRE 



"Super Burst" 



■ INTEGRAL Bright Bursts 

■ INTEGRAL Weak Bursts 
* RXTE Galloway et al. 2006 
□ SAS 3 Basinska et al. 1984 



10 

-2\ 



15 



E b (10 erg cm ) 



Fig. 6. Relation between peak flux, F pca k, and burst flu- 
ence, E\>. The dashed line indicates the lower limit of the 
observed photospheric radius expansion events. The two 
INTEGRAL bursts exhibiting PRE are labeled with the 
number 2 and 4 as reported in Tabled 



More than 200 bursts have been observed in 4U 1728- 

I I I j I 

34 since its discovery (Hoffman e t all 1197 6: White et al 
198fit iBasinska et allllQftdlFoster. Ross, fc Fabianl I l98f : 



Day fc Tawara 



Galloway et al 



.roster xvOl^ 
1 900J lFrancTl200lUvan StnmienL C : 
2003, 2006). Compared to the previously 



observed bursts, ours appear as ordinary bursts observed 
for this source. All the detected bursts have a rise time 
about 1 s with a e-fold decay time from 3-7 s in the 
3-20 keV energy band. Most of the bursts are observed 
up to 30-40 keV during the first 5-10 s, and the ener- 
getic ones are observed at the lowest inferred mass ac- 
cretion rate (see Fig. Such behaviour was observed 
e.g., also for GX 3+1, where a super burst and a unusu- 
ally long burst occurred at the lowest mass accretion r ate 
l(Kuulkers fc van der Krisl l2000UChenevez etaill200fift. 
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The averaged spectrum during the burst measured 
a blackbody temperatures from 1.8 to 2.4 keV with a 
bolometric flux from 1.7 x 10~ 8 to 4.5 x 10 -8 erg cm" 2 
s , respectively. This corresponds to an effective emis- 
sion radius (i.e., i? cff = di k P c F^ 2 T^ 2 cr" 1 / 2 (z + 1)~ 1/2 
where 2 + 1 — 1.31 is the gravitational redshift correction 
and a is the Stefan-Boltzman constant), i? e ffi from ~ 4.5 
d4.5kpc + 3.8 d4.5k pc km, assuming a NS radius of 10 km, a 
NS mass of 1.4 M Q and isotropic emission. This emission 
radius could be attributed to an extended part of the NS 
surface. 

In Fig. we show the bursts' peak flux, Fpcak, plot- 
ted as a function of the integrated burst flux, E\>, includ- 
ing the 32 and 104 bursts observed by SAS 3 and RXTE 
l|Basinska et allll984lGallowav et aTll200ri) . Note that a 
subset of the RXTE bu rsts were reported and discussed in 
jGallowav et all 12003^) . Most of the bursts with a higher 
peak flux as ~ 7 x 10~ 8 erg cm -2 s _1 show double-peaked 
behaviour, these energy-dependent double-peaked profile 
are also the result of PRE. All the bursts below the dashed 
line are bursts without PRE. The INTEGRAL bursts, ex- 
cluding bursts with PRE, occurred at a persistent emis- 
sion flux minimum and maximum level of ~ 2.2 x 10 -9 
and - 8.8 x 10~ 9 erg cm" 2 s" 1 (to - 2.4 x 10 3 ^ 9.4 x 10 3 
g cm' 2 s _1 ), respectively (3- 20 keV). The values w ith- 
out PRE observed by SAS 3 l|Basinska et all Il984|) are 
within the INTEGRAL observation. The averaged value 
of the persistent emission during the SAS 3 observation 
was ~ 3.25 x 10 -9 erg cm~ 2 s _1 (1-20 keV) and varied 
by a factor - 1.25. The INTEGRAL bursts with PRE, 
burst 2 and 4, occurred at the lowest mass accretion rate 
of ~ 1.7 x 10 3 g cm~ 2 s _1 . The bursts total fluence, and 
the bursts peak flux are anti-correlated with the mass ac- 
cretion rate. From a 30 days averaged RXTE/ ASM light 
curve (1.5-12 keV) llLevine et al.l El9961 the persistent flux 
of 4U 1728-34 shows a variation of about a factor three on 
a ~ 7.7 yr time scale and INTEGRAL observed 4U 1728- 
34 during a bright persistent emission phase (see Fig. EJ . 
Indeed the only two bursts with PRE where observed at 
the the lowest mass accretion rate at the start of the 
INTEGRAL observation, all other bursts without PRE 
are observed at higher mass accretion rate. 

The lowest inferred mass accretion rate at the ~ 7 x 
10~ 8 erg cm~ 2 s _1 line (see Fig. is m « 2.4 x 10 3 g 
cm~ 2 s _1 . From the observed burst properties (see Table 
^) and inferred mass accretion rates, the present theory 
predicts that these bursts ar e pure helium burning (e.g. 
IStrohmaver fc Bildstenl . r2004|) . For helium flashes, the fuel 
burns rapidly, since there are no slow weak interactions, 
and the local Eddington limit is often exceeded. These 
conditions lead to PRE bursts with a duration, set mostly 
by the time it takes the heat to escape, of the order 
of 5-10 s, as observed. The INTEGRA L-observed bursts 
without PRE with a mass accretion rate of > 2.4 x 10 3 
g cm~ 2 s _1 could also arise from helium burning. In 
the framework of the thermonuclea r-flash models (e.g., 
iLewin. van Paradiis. fc Taaml Il995|) the burst duration, 
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Fig. 7. 30-day averaged RXTE/ ASM light curve in the 
1.5-12 keV energy band showing the long term variability 
of 4U 1728-34 from 1996 to 2005. The vertical dashed lines 
indicate the time interval of the INTEGRAL observations. 
The ASM count rates has been converted into flux using 
1 Crab Unit for 75 cts/s. The arrows indicate the time of 
the type I bursts. 



t = Sb/Fp ea k < 10 s, and the ratio of the observed per- 
sistent flux to the net peak flux, 7 = F pers / Fp ea ki indi- 
cate a hydrogen-poor burst, i.e., pure helium bursts. We 
found a good agreement with the predicted amount of lib- 
erated helium fuel during the bursts with the observed 
accreted mass between two bursts (calculated only for the 
bursts with known a parameter). The amount of fuel liber- 
ated in the thermonuclear burning is calculated as Lb/eHe, 
where Lb is the observed burst total fluence and ehc ~ 1.7 
Mev/nucleon « 1.6 x 10 18 erg g _1 is the He energy re- 
lease. For instance, for burst 2, AM = MAt ~ 1.2 x 10 21 
g is in agreement with the predicted amount of liber- 
ated helium fuel during this burst Lb/enc ~ 1-3 x 10 21 
g. We observed pure helium bursts during the different 
spectral states. 4U 1728-34 shows the same burst prop- 
erties as the ultra compact 10 minute binary 4U 1820-30 
l|Cummind . 120031). In Fig. El we also s how the so-called 
"super burst" bv iBasinska et alJ J^84). Super bursts, as 
currently understood and observed, have a duration of 
several hours, a total fluence of ~ 10 42 erg and are be- 
lieve d to be due to unstable burning of a thick carbon layer 
(e.g. IStrohmaver fc Bildstenl 120 04). The observed "super 
burst" bv iBasinska^T^lJ l|l984|) was most likely also a 
helium burst with PRE triggered after a longer burst re- 
currence time. 



5. Summary 

We analyzed the spectral and bursting behaviour from 
the simultaneous INTEGRAL IBIS/ISGRI and JEM-X 
observation of 4U 1728-34, which was spatially well distin- 
guished from the neighbouring source, the Rapid Burster. 
The broad band spectra from 3-200 keV is well described 
by a thermal Comptonization model with seed photons 
from the neutron star surface scattered in a shock-heated 
accretion column above the hot spot region, or from an op- 
tically thick boundary layer between accretion disc and NS 
surface. The thermal blackbody soft emission could arise 
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from the accretion disc. The spectral evolution observed 
from an intermediate/hard state to the soft state are en- 
tirely typical of atoll sources. With increasing mass accre- 
tion rate the source spectra become brighter and softer. 
We observed a decreasing electron temperature from 35-3 
keV an increase of the optical depth from 0.5-5, during 
the source luminosity increase from 2-12% of Eddington, 
i.e. from the intermediate/hard state to the soft state. 

The 36 type I bursts show an anti-correlation between 
burst peak flux or fluence with the mass accretion rate. 
At low mass accretion rate, intermediate/hard state, most 
of the bursts show PRE with a peak flux exceeding the 
Eddington luminosity, at higher mass accretion rate, soft 
state, the bursts were normal type I bursts. We observed 
normal X-ray bursts for this source involving pure helium 
cither during PRE in the intermediate/hard state, or dur- 
ing the soft state. 
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